How big is an outbreak likely to be? Methods for epidemic final-size calculation

Epidemic models have become a routinely used tool to inform policy on infectious disease. A particular interest at the moment is the use of computationally intensive inference to parametrize these models. In this context, numerical efficiency is critically important. We consider methods for evaluating the probability mass function of the total number of infections over the course of a stochastic epidemic, with a focus on homogeneous finite populations, but also considering heterogeneous and large populations. Relevant methods are reviewed critically, with existing and novel extensions also presented. We provide code in Matlab and a systematic comparison of numerical efficiency.Thomas House, Joshua V. Ross and David Sir

A discrete MMAP for analysing the behaviour of a multi-state complex dynamic system subject to multiple events.

A complex multi-state system subject to different types of failures, repairable and/or nonrepairable, external shocks and preventive maintenance is modelled by considering a discrete Markovian arrival process with marked arrivals (D-MMAP). The internal performance of the system is composed of several degradation states partitioned into minor and major damage states according to the risk of failure. Random external events can produce failures throughout the system. If an external shock occurs, there may be an aggravation of the internal degradation, cumulative external damage or extreme external failure. The internal performance and the cumulative external damage are observed by random inspection. If major degradation is observed, the unit goes to the repair facility for preventive maintenance. If a repairable failure occurs then the system goes to corrective repair with different time distributions depending on the failure state. Time distributions for corrective repair and preventive maintenance depend on the failure state. Rewards and costs depending on the state at which the device failed or was inspected are introduced. The system is modelled and several measures of interest are built into transient and stationary regimes. A preventive maintenance policy is shown to determine the effectiveness of preventive maintenance and the optimum state of internal and cumulative external damage at which preventive maintenance should be taken into account. A numerical example is presented, revealing the efficacy of the model. Correlations between the numbers of different events over time and in non-overlapping intervals are calculated. The results are expressed in algorithmic-matrix form and are implemented computationally with Matlab.Junta de Andalucía, Spain, under the grant FQM307Ministerio de Economía y Competitividad, España, MTM2017-88708-PEuropean Regional Development Fund (ERDF

Optimizing Performance of Continuous-Time Stochastic Systems using Timeout Synthesis

We consider parametric version of fixed-delay continuous-time Markov chains (or equivalently deterministic and stochastic Petri nets, DSPN) where fixed-delay transitions are specified by parameters, rather than concrete values. Our goal is to synthesize values of these parameters that, for a given cost function, minimise expected total cost incurred before reaching a given set of target states. We show that under mild assumptions, optimal values of parameters can be effectively approximated using translation to a Markov decision process (MDP) whose actions correspond to discretized values of these parameters

Distributed Synthesis in Continuous Time

We introduce a formalism modelling communication of distributed agents strictly in continuous-time. Within this framework, we study the problem of synthesising local strategies for individual agents such that a specified set of goal states is reached, or reached with at least a given probability. The flow of time is modelled explicitly based on continuous-time randomness, with two natural implications: First, the non-determinism stemming from interleaving disappears. Second, when we restrict to a subclass of non-urgent models, the quantitative value problem for two players can be solved in EXPTIME. Indeed, the explicit continuous time enables players to communicate their states by delaying synchronisation (which is unrestricted for non-urgent models). In general, the problems are undecidable already for two players in the quantitative case and three players in the qualitative case. The qualitative undecidability is shown by a reduction to decentralized POMDPs for which we provide the strongest (and rather surprising) undecidability result so far

Petri Nets Validation of Markovian Models of Emergency Department Arrivals

International audienceModeling of hospital’s Emergency Departments (ED) is vital for optimisation of health services offered to patients that shows up at an ED requiring treatments with different level of emergency. In this paper we present a modeling study whose contribution is twofold: first, based on a dataset relative to the ED of an Italian hospital, we derive different kinds of Markovian models capable to reproduce, at different extents, the statistical character of dataset arrivals; second, we validate the derived arrivals model by interfacing it with a Petri net model of the services an ED patient undergoes. The empirical assessment of a few key performance indicators allowed us to validate some of the derived arrival process model, thus confirming that they can be used for predicting the performance of an ED

Delayed Capital Injections for a Risk Process with Markovian Arrivals

In this paper we propose a generalisation to the Markov Arrival Process (MAP) risk model, by allowing for a delayed receipt of required capital injections whenever the surplus of an insurance firm is negative. Delayed capital injections often appear in practice due to the time taken for administrative and processing purposes of the funds from a third party or the shareholders of a firm. We introduce a MAP risk model that allows for capital injections to be received instantaneously, or with a random delay, depending on the amount of deficit experienced by the firm. For this model, we derive a system of Fredholm integral equations of the second kind for the Gerber-Shiu function and obtain an explicit expression (in matrix form) in terms of the Gerber-Shiu function of the MAP risk model without capital injections. In addition, we show that the expected discounted accumulated capital injections and the expected discounted overall time in red, up to the time of ruin, satisfy a similar integral equation, which can also be solved explicitly. Finally, to illustrate the applicability of our results, numerical examples are given

A general piecewise multi-state survival model: Application to breast cancer

Multi-state models are considered in the field of survival analysis for modelling illnesses that evolve through several stages over time. Multi-state models can be developed by applying several techniques, such as non-parametric, semi-parametric and stochastic processes, particularly Markov processes. When the development of an illness is being analysed, its progression is tracked periodically. Medical reviews take place at discrete times, and a panel data analysis can be formed. In this paper, a discrete-time piecewise non-homogeneous Markov process is constructed for modelling and analysing a multi-state illness with a general number of states. The model is built, and relevant measures, such as survival function, transition probabilities, mean total times spent in a group of states and the conditional probability of state change, are determined. A likelihood function is built to estimate the parameters and the general number of cut-points included in the model. Time-dependent covariates are introduced, the results are obtained in a matrix algebraic form and the algorithms are shown. The model is applied to analyse the behaviour of breast cancer. A study of the relapse and survival times of 300 breast cancer patients who have undergone mastectomy is developed. The results of this paper are implemented computationally with MATLAB and R.Ministerio de Economía y Competitividad FQM-307European Regional Development Fund (ERDF) MTM2017-88708-PUniversity of Milano-Bicocca 2014-ATE-022